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Abstract 

A generalized inverse seesaw model, in which the 9x9 neutrino mass matrix has vanishing 
(1,1) and (1,3) submatrices, is proposed. This is similar to the universal two-zero texture which 
gives vanishing (1,1) and (1,3) elements of the 3x3 mass matrices in both the charged lepton and 
neutrino sectors. We consider the Z G x Z & group to realize such texture zeros in the framework 
of the generalized inverse seesaw model. We also analyze the universal two-zero texture in the 
general case and propose two ansatze to reduce the number of free parameters. Taking account 
of the new result of 9 13 from the Daya Bay experiment, we constrain the parameter space of 
the universal two-zero texture in the general case and in the two ansatze, respectively. We find 
that one of the ansatze works well. 
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1 Introduction 



The canonical seesaw mechanism [T] is successful in generating small masses of left-handed neutri- 
nos, but it has no direct experimental testability and encounters a potential hierarchy problem [2]. 
In the type-I seesaw model with heavy right-handed neutrinos iV R , the left-handed neutrinos v L 
can gain small masses M v ~ M-^M^M^ thanks to the huge right-handed neutrino masses M R . 
However, to obtain M v ~ O(0. 1) eV, one has to require M R ~ 0(1O 14 ) GeV, if M D is assumed to 
be at the electroweak scale (~ O(10 2 ) GeV). This makes the right-handed neutrinos far beyond 
the detectability of any colliders. The hierarchy problem is that a very high seesaw scale will lead 
to large corrections to the Higgs mass, which makes the Higgs mass of the order of the electroweak 
scale unnatural. The inverse seesaw model |3j can solve these problems. Moreover, it is possible to 
predict light sterile neutrinos naturally |4J and provide rich phenomenology such as the non-unitary 
effect and leptogenesis [5]. 

The generalized inverse seesaw model (GISM) is an extension of the canonical seesaw mech- 
anism by introducing three right-handed neutrinos iV Ri (for i = 1,2,3), three additional gauge- 
singlet neutrinos 5 Ri and a scalar $ into the standard model (SM). The Lagrangian in the charged 
lepton and neutrino sectors [6] is written as 

-A = ^^R + ^D^R + ^S^R + ^M R 7V R+ l^M^ R + h.C., (1) 

in which H, £ L and E R stand for the Higgs doublet, three lepton doublets and three charged-lepton 
singlets, respectively in the SM and H = ia 2 H*. Here Y^ and Yg are 3x3 Yukawa coupling 
matrices, and M R and are 3x3 symmetric Majorana mass matrices. After spontaneous 
symmetry breaking (SSB), the scalars acquire their vacuum expectation values (VEVs), and we 
gain the 3x3 charged lepton mass matrix M l = Y,v(H) / \/2 and the 9x9 neutrino mass matrix 
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in the flavor basis, in which M D = Y D v(H)/y/2 and M s = Y s v(§)/yj2. Here v(H) and f($) are 
the VEVs of H and respectively. The GISM degrades to the original inverse seesaw model 
(OISM) when M R = is taken. It can also accommodate a larger range of the sterile neutrino 
masses than the OISM [3]. 

If we regard each submatrix of Ai in Eq. (2) as a complex number, we turn to a typical pattern 
of two-zero textures Different from the models given in Ref. [7], where M ; is chosen to be 
diagonal and only M v has the two-zero texture, we propose the universal two-zero texture (UTZT) 
[8], in which both M l and M v have two-zero textures. As the similar texture zeros of quark mass 
matrices can interpret the smallness of flavor mixing angles in the quark sectors [9], we expect the 
UTZT will give us a better understanding of the lepton flavor mixing. We write out the charged 
lepton and left-handed neutrino mass matrices as 
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Some work of this texture has been done in Ref. [10\ [TTj . Generally, texture zeros can be obtained 
from Abelian symmetries |12j . Later we will see that by means of these symmetries the two-zero 
texture of can be directly derived. For the light left-handed neutrino matrix M u , if M D , M R , 
Mg and M all have the two-zero textures, which will be a natural result from the symmetries, 
the seesaw mechanism can guarantee that M v achieves the two- zero texture [H [11] . 

Recently, the Daya Bay collaboration reported a relatively large # 13 [13] with its best-fit (±1<t 
range) value # 13 ~ 8.8° ± 0.8°. It is confirmed by the RENO experiment [2]. The experimental 
results of large # 13 give us two motivations for the UTZT. (1) Two phenomenological strategies 
towards understanding lepton flavor mixing are outlined in Ref. [15] : the first one is to start from 
a nearly constant flavor mixing pattern, and the second one is to associate the mixing angles with 
the lepton mass ratios. While it is a nontrivial job to generate a large 9 13 from the first strategy 
according to flavor symmetries, one may pay more attention to the second strategy. To implement 
the second strategy, one generally requires some elements of M l and M v to be zeros or sufficiently 
small compared with their neighbors, and the two-zero texture is a typical example of this kind. 
(2) As discussed in Ref. [16], where M l is diagonal and M v has a two-zero texture, it is more likely 
to obtain a large # 13 if M v has texture zeros as in Eq. (3) compared with the other texture zeros. 
Taking advantage of this kind of texture zeros, we expand our discussion to the scenario that both 

and M u have such texture zeros. We expect that such texture can also gain a large # 13 . 

The rest of this paper is organized as follows. In section 2, we propose a model to connect 
the GISM with the UTZT under the discrete Abelian group Z 6 x Z 6 . With this model, we can 
realize the two-zero textures of M D , M R , M g and M^. However, the realization of the two-zero 
texture of M u is a little non-trivial. Section 3 is devoted to see how the two-zero texture of M v is 
realized. In section 4, the UTZT is used to explain the lepton flavor mixing, especially for large 
9 13 . Both analytical and numerical results are presented. The predictions for the effective masses 
in the tritium beta decay and neutrinoless double beta (Oz/2/3) decay are also given in this section. 
Since the UTZT in the general case has several adjustable parameters, it does not get stringent 
experimental constraints. In section 5, we consider two ansiitze of the UTZT to constrain the 
parameter space. Ansatz (A) is a natural approximation based on our model built in section 2, 
and ansatz (B) is a special case which has been considered in Ref. [8]. Section 6 is the conclusion 
of our paper. 

2 A model connecting the UTZT with the GISM 

In this section, we illustrate a way to connect the GISM with the UTZT. We rewrite the Lagrangian 
in the charged lepton and neutrino sectors as 

+\Wi(YR^X a N K] + \Wi(y^4> a S K] + h.c. , (4) 

in which the repeated indices are summed. In our model, we introduce three scalars into each 
term, so a = 1,2,3. Comparing Eq. (4) with Eq. (1), we can see that some replacements have 
been done. Yfl, Y B H, Y s & are replaced by YfH -, YgH a , Y g a <I> a , respectively, and the scalars 
X a , 4> a are introduced to give the Majorana masses of N R , S R , respectively. The purpose to do 
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Table 1: The charges of the fermions and scalars under Z f 







E Ri 










x a r 







4 


2 


1 


4 


3 


2 4 



Table 2: The charges of the fermions and scalars under Z ( 
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these replacements has nothing to do with the GISM but to give the two-zero textures of the mass 
matrices M l} M D , M s , M R and M^. 

A model for connecting the GISM with the UTZT can be built based on a direct product of 
groups G lQi x G 2q2 = Q: 

• Each fermion or scalar transforms under the group G 1 with a charge q 1 . This rule aims to 
realize the GISM. Since it is flavor-blind, different flavors in the same multiplet (e.g., N { and 
Nj with i / j) have the same charges q x , and different scalars in the same Yukawa coupling 
(e.g., H a and H b with a^b) have the same charges q l , too. 

• Each fermion or scalar transforms under group G 2 with a charge q 2 - We choose G 2 to be 
an Abelian group Z n to give the UTZT. In this case, different flavors in the same multiplet 
should have different charges q 2 , and different scalars in the same Yukawa coupling term 
should also have different charges q 2 . 

Generally speaking, there are many possibilities to choose G l and G 2 , and it is essentially unnec- 
essary to require that they be equal to each other. Nevertheless, in view of the similar structures 
of M and M x u , we assume G l = G 2 = Z n . 

In our model, we choose n = 6 and Q = Z eq x Z eq . The discrete Abelian group Z 6 is given 
by Zg = {1, ui , uj 2 , uj 3 , w 4 , u} 5 } , where uj = e" 1 "/ 3 . In Tables 1 and 2, we list the charges q 1 and q 2 
for each field, respectively. The invariance of the Lagrangian under the Z 6 x Z 6(?2 leads to the 
following textures of the Yukawa coupling matrices: 

/ \ 

and F| ~ x , (5) 
\ j 

for Y% = Y t a , Yg, Y s a , Y^ and Y" M a . After SSB, the scalars gain their VEVs, and we are left with 
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the mass terms 
-L, 



£ L M^ R + ^M D 7V R + N§,M S S R + ^M R N R + ^ R M^ R + h.c. 
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where Mj, M D , M g , M R and M^ are mass matrices originating from the Yukawa coupling matrices 
and VEVs of the scalars. Taking M D for example, we arrive at 



Mi 
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-j= [yMh 1 ) + yMh 2 ) + yMh 3 )} , 



(7) 



in which v(H a ) is the VEV of H a . All the mass matrices Mj, M D , Mg, M R and M^ have the same 
two-zero textures as 

x \ 
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In appendix A, we show that the mass matrix of light left-handed neutrinos is given by a seesaw-like 
formula in the physical region: 



M v = -M D (M R -M s M~ 1 Mi) l M?. 



(9) 



With this formula, one can prove that M v also follows the two-zero texture as in Eq. (8) |17t [TT] . 
A detailed analysis will be given in the next section. 

We remark that besides Z 6 , lots of discrete Abelian groups Z n can connect the GISM with the 
UTZT. Even under the same discrete Abelian group, a different arrangement of the charge q 2 may 
cause different textures of the Yukawa coupling matrices Y\, Y\ and Y|, but it keeps the textures 
of mass matrices as in Eq. (8) unchanged. In brief, there are many possibilities to link the GISM 
with the UTZT. However, if one requires that the Abalian discrete symmetry be anomaly-free, 
one must pay attention to the arrangement for the charges q l and q 2 of each field to guarantee the 
anomaly- free conditions [18J. Then some arrangements for the charges q 1 and q 2 will be ruled out. 



3 The mass texture of active neutrinos 

We have proposed a way to realize the two-zero textures of Mj, M D , M R , M s and M„. These 
textures can be obtained immediately from flavor symmetries under the direct product of discrete 
Abelian groups. However, a realization of the two-zero texture of M v is not so obvious. To find its 
texture, we must turn to the matrices M D , M R , M s and M„, all of which have the same texture 
zeros. In a way similar to the proof in Refs. |17| and [TT], after giving the two-zero textures of 
M D , M R , Mg and M„, we can prove that the two-zero textures manifest themselves again in M u , 
as a consequence of Eq. (9) . 

We express each matrix M a (for a = D, S, //) as 
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It is easy to find the inverse matrix of has another type of texture zeros 
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Then, using the seesaw formula M x 



-MgM 1 M§ , we find M x has the two-zero texture as 
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Thus the two-zero texture is invariant under the seesaw transformation. 

Repeating the above process for M v = —M D (M R + M x )~ l M^, we finally obtain that M v has 
the two-zero texture as in Eq. (3). The non-zero entries are given by 
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It is an exact consequence of the GISM and two-zero textures of M D , M R , M s and M„. 

All the 3x3 mass matrices Mj, M u , M D , M R , M s and has parallel structures with each 
other. And they are all fractally similar to the 9x9 GISM neutrino matrix A4. These similarities 
can be guaranteed in the framework of flavor symmetries. 



4 Flavor mixing in the UTZT 

In this section we analyze the flavor mixing in the general UTZT case. The renormalization-group 
effect might in general modify the two-zero textures of M l and M v , but it is negligibly small in the 
inverse seesaw model [19] because the TeV seesaw scale is so close to the electroweak scale. Hence 
we just discuss the UTZT at the electroweak scale. 
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The charged lepton and left-handed neutrino mass matrices with two-zero textures have been 
given in Eq. (3), where A t v , B [ v , C l v and D l v are complex numbers. Some works on this texture 
have been done in Refs. [8] and |10j . but a general analysis has been lacking in the literature. 

As a symmetric matrix, M t can be diagonalized as M ; = V l M l Vj r . Here M t = Diag{m e , m , m T }, 



V l = Qptf, Q l = Diag{e' 
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jPi, l},P t = Diag{e^ 
( c 



i 









v- 



e l7r } and U l is given by 
ft 



\ 







(15) 



in which c a = cos9 a , s a = sin 6^ (for a = e,fj,,r) and 

Similarly, M v can be diagonalized as M v = V u M v Vj . Here M v 



Diag{m l5 m 2 , m 3 }, V u 



Q V U V P V , Q v = Diag{e* 
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The Maki-Nakagawa-Sakata-Pontecorvo (MNSP) matrix [20J is defined byV = V{V l 



in which Q = Diag{e* a , , 1} and a, j3 are two combined parameters defined as a 
/3 = j3 v — /3j, respectively. V can be parametrized as V = QUP. Here 

/l \ / c™ st,\ / d 



(16) 



p}u}qu u p v , 



at, 



U 







-23 



, 23 



'23 c 23/ 



v- 



-13 



" s 13 



12 



'12 



12 



12 



(17) 



V 



in which Cj ■ = cos 4 ■, = sinO- (for = 12, 23, 13) and s 13 = s 13 e* <5 . P and Q are two diagonal 
phase matrices. As the charged leptons are the Dirac fermions, Q is unphysical and can be rotated 
away by the phase redefinition of the charged lepton fields. But for the Majorana neutrinos, only 
one overall phase in P can be rotated away and the other two phases are physical. In this case, P 
can be parametrized as P = Diag{e ip , e lCT , 1}. 

4.1 Charged leptons 

Here we derive some relations of the mixing parameters in the charged lepton sector. Since the 
(1,1) and (1,3) elements of M h are equal to zeros, we obtain 
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A straightforward calculation leads us to the relations of the angles 

cot 2 9 e = s 2 (J xfyf cot 4 Ofj, - sin 2 5^ - cos <0 \ J V? cot 4 0^ - sin 2 ^ + cos SA , 
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and those of the phases 
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7^-7. + V 2 = °. ( 2 °) 

where x ; = m e /m„ and j/j = m fl /m T . Taking m e = 0.486 MeV, = 102.7 MeV and m T = 
1746 MeV at the electroweak scale [21] as inputs, we get x h = 0.0047 and y l = 0.059. To assure 
that Eq. (19) have a real and positive solution, we require 

0^9^^ arctan y/x^ ~ 1° , 
9 T < arctan rj 6° , 

< 9 e ^ 90° . (21) 

In particular, T « arctan^ w 4° for <^ = ±90°, < 9 T < 4° for \S^\ < 90°, and 4° < T 6° 
for |<5J ^ 90°. Due to the large mass hierarchy of the charged leptons, 0^ and 9 T are very small. 
They can be regarded as the corrections to the MNSP matrix. Suppressed by s„, the phase 5^ has 
little influence in the MNSP matrix. Particularly, we have three special cases: 
(1) tan0 e < l/y/V v 

taxi « y/xi y l tan9 e , 
tan 9 T « , 

7e «7 M ±90°. (22) 



(2) tan0 e > 1/^, 



tan » 



tan 9 T «s . / — cot e , 
V 

7 e «7 M - V2±90°. (23) 

(3) tan0 e ~ (l/y/yTJ, one can find tan0 M ~ O {y/x~jy[) and tan0 r ~ O {y/xj) from Eq. (19). 
In the leading-order approximation of and s T , we obtain 

2 _ -x^fcos(7 e -7 M ) 
S/1 ~ x l + y { tan 2 9 e 
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-x l cos(7 e - 7 M ) 
X; + y ; tan 2 9 e 



sin 5 M ss x lVl sin(0 e - 6^) , (24) 

and 7 e — 7^ is arbitrary. 
4.2 Neutrinos 

For the left-handed neutrinos, since the (1,1) and (1,3) elements of M v equal zeros, we obtain 
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Later in the numerical calculations, we will see that 9 2 is a small angle, in the same order of 
magnitude as 9 l3 . In this case, we find m 1 < m 3 from Eq. (25). Only the normal hierarchy of 
neutrino masses is possible in this texture. A straightforward calculation leads us to the relations 
of the angles 

cot 2 9 1 = s\ (yj x 2 v yl cot 4 9 2 — sin 2 5 2 — cos 5 2 ^j {^\J vl c °t 4 $2 ~~ s ™ 2 ^2 + cos ^ 



tan 2 # 3 = 



yj xlyl cot 4 9 2 — sin 2 5 2 — cos S 2 



2 cot 4 9 — sin 2 <5 9 + cos S 



(26) 
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and those of the phases 
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(27) 



in which x v = m 1 /m 2 , y v = m 2 /m 3 and x u , y u < 1. To make Eq. (26) have a real and positive 
solution, we require 



^ 9 2 ^ arctan yjx v y v , 
^ 9 3 ^ arctan y/2x l/ , 
s$ ^ s: 90° . 



(28) 



4.3 The MNSP matrix 

We have obtained some relations of the mixing parameters in both the charged lepton and left- 
handed neutrino sectors. Using these parameters, we can calculate the mixing angles in the MNSP 
matrix and some other physical observables. And using the experimental constraints, we may find 
the allowed ranges of the parameters and make predictions for the observables. 

The MNSP matrix V can be calculated through V = V l V v . Considering the smallness of s^, 
s T and s 2 , we obtain the approximate expressions of the mixing angles 9 13 , 9 12 and 9 23 : 
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These expressions hold to the first order in s^, s T and s 2 . We make some comments on the formulas 
of the mixing angles in Eq. (29): 



• Note that 9 13 is in the same order of magnitude as 9 2 , and 9 2 ^ arctan ^Jx v y v = arctan y / m 1 /m 3 . 
To generate a relatively large 9 13 , m 1 cannot be too small. 

• Since 9^ and 9 T are small, # 12 ~ 9 3 holds. The two-zero texture in the charged lepton sector 
just has a small contribution to 9 12 . 
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• # 2 3 is an overall result of 9 1 , 9 e and j3. The two-zero textures in both the charged lepton and 
neutrino sectors may have large contributions to 9 23 . 

We conclude that except for 6> 12 , both 9 13 and 8 23 may receive relatively large corrections from the 
charged lepton sector. This is one of the features that make the UTZT different from the texture 
zeros discussed in Ref. [7], in which M l is diagonal and only M u has texture zeros. 

The strength of CP violation in the neutrino oscillation experiments is measured by the Jarlskog 
invariant J = In^V^V^V^^Ti) = c i2 s i2 c 23 s 23 c i3 s i3 sm ^ E2]- For current experimental data of 
8 13 , one may expect a relatively large J if the CP-violating phase b is not suppressed. In the 
leading-order approximation of s„, s T and s 2 , 



J ~ s J + s T J T + s 2 J 2 , (30) 



in which 



J M = - S 3 C 3 [ S l C e sin (« + <*/J - C l S e sin ( a + S fJ ,~ P)] {44 + 44 + 2 C 1 S 1 C £ S e COS (3) , 

J T = — S3C3 [s 1 s e sin a + c 1 c e sin (a — /?)] {44 + s\4 — 2c 1 s 1 c e s e cos 0) , 

^2 = s 3 c 3 (c e s e sin /3 cos <5 2 + c 1 s 1 sin 5 2 cos 26* e — c e s e cos /3 sin <5 2 cos 20]^) . (31) 

One can see that the first term s^J^ is in general the smallest one because of the smallness of s^, 
and the last two terms s T J T and s 2 J 2 may have comparable contributions to J . 

The 0^2/3 decay experiments are important for examining if neutrinos are the Majorana 
fermions. One key parameter in such experiments is the effective mass (m) ee = (VM U V T ) 11 = 
(V^ M u V l *) 11 . The pattern in which M v has the two-zero texture in Eq. (8) and M l is diagonal 
gives {m) ee = 0. Different from such a pattern, the UTZT that we are considering here yields a 
non-zero {m) ee . In the leading-order approximation of s^, s T and s 2 , {m) ee reads 

(m) ee « 2|(^)* 1 (MJ 21 | « 2\c e s T + 3,8^1, (32) 

in which 

\A V \ « Imgs^ae 2 ^ - m 1 c 1 c 3 s 3 e 2171 + maCiCgSge 2 * 72 ! . (33) 
In comparison, the effective mass (m) e in the tritium beta decay is given by 



(m) e = yJ(VM^) u « L4J. (34) 

Then, we arrive at 

{m) " « 2|c eSr + Se ^|. (35) 



Me 

If we assume # T = 4° and ignore the smallness of 9^, then we obtain {'m) ee /(m) e ~ 0.1. 
4.4 Numerical results 

In the numerical calculations, we choose 7 free parameters 9 e , 9 2 , <L, <5 2 , a, /3 and m 1 as inputs. 
The values of the charged lepton masses have been given in section 4.2. To be compatible with the 
experimental results, we choose Am 21 ~ (7.4 — 7.8) x 10~ 5 eV 2 , Ara^ ~ (2.4 — 2.7) x 10~ 3 eV 2 , 
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<9 23 ~ (42° - 49°), 9 12 ~ (33° - 35°) and the new data 13 ~ (8.0° - 9.6°) from the Daya Bay 
experiment as constraints. With the help of these data, we can obtain the allowed ranges of the 
input parameters and calculate the observables. 

In Fig. 1, we show the comparison of the values between 9 1 and 9 e and that of the values 
between 9 2 and 9 T . The first two angles are dominant parameters in the expression of 9 23 , while 
the last two angles are dominant parameters in the expression of 13 (see Eq. (29)). Numerically, 
we obtain 9 1 ~ (24° - 72°) versus e ~ (0 - 90°) for 23 , and 9 2 ~ (4° - 13°) versus 9 T ~ (0 - 6°) 
for 9 13 . A lot of points are located around 9 T = 4°, indicating that <L « ±90° is favored. 

In Fig. 2, we show the parameter space and some phenomenological predictions in the general 
case. We plot the allowed regions of (0 13 ,m 1 ) and {9i 2 ,9 23 ) parameters first in the figure. The 
lightest neutrino mass m 1 is constrained in the range (0.001 — 0.015) eV. The points of the mixing 
angles 9 12 , 9 23 and 9 13 are nearly evenly distributed in the full parameter space. Predictions for 
parameters related to CP violation are shown then. There is little restriction on the combined 
input parameters a and f3 except that /3 is more likely to approach ±90°. For the Majorana phases 
p and a, the relation p ~ a ± 90° holds roughly. The numerical result of the Jarlskog invariant J 
is also shown in Fig. 2. Due to the largeness of 9 l3 , \ J\ can reach several percent. Concretely, it 
can maximally reach 0.03 at 13 = 8° and 0.04 at 13 = 9.6°. The effective masses in the tritium 
beta decay and 0v2(3 decay are shown at the end of Fig. 2. One can see that the ratio (m) ee /{m) e 
is of 0(0.1) in most cases. Since {m) e ~ 0.01 eV is referred in Fig. 2, (m) ee can maximally 
reach 10 -3 eV. However, this is still below the sensitivity of the near future experiments, which is 
expected to be (m) ee ~ (1 - 5) x 10~ 2 eV [23]. 

One can reconstruct the charged lepton and left-handed neutrino mass matrices with the help 
of the experimental constraints. Considering that there are cancelations in some special cases, 
leading to vanishing values of A lv , B lv , C lv or D l v , the positive lower bounds may not exist. 
But one can expect that there are some ranges in which most of the points are located. In our 
calculation, we find that 95% of the points are located in the following ranges: 

\A t \ ~ (7.4 - 31) MeV , \B t \ ~ (0.046 - 0.94) GeV , 

|C,| ~ (0.96 - 1.8) GeV , \D t \ ~ (0.11 - 1.8) GeV , (36) 

and 

\A„\ ~ (0.0073 - 0.018) eV , \B V \ ~ (0.019 - 0.028) eV , 

\C V \ ~ (0.011 - 0.040) eV , \D V \ ~ (0.010 - 0.048) eV . (37) 

In the neutrino sector, all the elements of M u are in the 0(0.01) eV order. But in the charged 
lepton sector, the elements of M ; vary within some wide ranges because of the uncertainty of 9 e . 



5 Large 13 and two ansatze of the UTZT 

In the previous section, we have considered the UTZT in the general case. Since there are 7 free 
parameters as inputs, it does not get stringent experimental constraints. We shall consider some 
special cases of the UTZT to simplify its texture. 
First, we assume that the condition [8 J 

arg(C^) + arg(ZV) = 2 arg(i^) (38) 



11 





( ° 


A,J 
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1*1, J 


\D l>u 



is satisfied. Then M l and M v are respectively decomposed into 

M^PjMf^r and M v = P?M v P v e 2 ^ , (39) 

in which 

/ \\ u \ \ 

(40) 

In the following discussions, we turn to two different ansatze: ansatz (A), \A lu \ = \D lu \; and 
ansatz (B), \C t \ = \B h \ and \C V \ = \D V \. 

5.1 Ansatz (A) 

We propose to consider this new ansatz, in which both \Aj\ = \Dj\ and \A V \ = \D V \ hold. Our 
motivations are based on the model which we built in section 2: 

• In the charged lepton sector, the (1,2) and (3,3) entries of the Yukawa coupling matrix YJ 1 
are nonzero. It is natural to assume that they have the same magnitude: | (l^ 1 ) 12 1 = | (Y^ 1 )33 1 - 
After SSB, we arrive at | ( A^ ) 12 1 — | C ) 33 1 , or equivalently, |^4J — \Dj\- In this assumption, 
we can reduce the number of free input parameters. This equality can be realized in the 
non-Abelian discrete group A 5 [Mj with suitable arrangements of the particle content^}]. 

• Applying the above discussion to the neutrino sector, we are led to 

|(M D ) 12 | = |(M D ) 33 | , |(M S ) 12 | = |(M S ) 33 | , 

|(M R ) 12 | = |(M R ) 33 |, |(M M ) 12 | = |(M M ) 33 | . (41) 



Then, using the inverse seesaw formula in Eq. (9), we arrive at 

|(M D ) 12 | 2 |(M M ) 12 | 
|(M s ) 12 |2-|(M m ) 12 ||(M r;i2 



\(M V ) 12 \ = \(M V ) 33 \ = u ^) {M ? ]l f} M ^l ; , , (42) 



or equivalently, \A V \ = \D U \. 

In ansatz (A), the mass matrices in both the charged lepton and left-handed neutrino sectors 
can be solved exactly in terms of their mass eigenvalues. In the left-handed neutrino sector, we 
have the expression of M v in terms of its three mass eigenvalues 

\A V \ = (m 1 m 2 m 3 ) 1/3 , 

\B U \ = [(m 1 ?7i 2 m 3 ) 1 ^ 3 (m 1 — m 2 + m 3 ) — 2 (m 1 m 2 m 3 ) 2 ^ 3 

+m l m 2 — m 1 m 3 + m 2 m 3 ] 1 ^ 2 , 
\C V \ = m 1 — m 2 + m 3 — (m 1 m 2 m 3 ) 1 ^ 3 , (43) 



t We may arrange i L and E R as the triplets, H 1 as a singlet, and embed H 2 and H 3 to a 5-plet in A 5 . After 
H 1 gains its vacuum expectation value, we are led to A l = D l . With a suitable vacuum alignment for the 5-plet, no 
additional mass term will be introduced and the two-zero texture is preserved. 
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and that of U u in terms of the ratios of the eigenvalues 

x uVu Q'u)Q'v 

U u = k ul (x u y u - a v )x v y v -k u2 (y u + a v )y u k u3 (l - a u ) 

\ k ul b v x u y u k v2 b v y v k v3 b v J 

where a v = \A u \/m 3 , b u = \B v \/m 3 , c v = \C u \/m 3 and 

Ki = [{ a l + x lvl) ( x uVu ~ a v? + xly 2 v bl] " 1/2 , 



(44) 



V2 



k 



;/3 



[(al + yl){y u + af +ylb 
[(al + l)(l-a u f + bl] 



2,21-1/2 



-1/2 



(45) 



In the charged lepton sector, after replacing the index v — > I and the masses (m 1 ,m 2 ,m 3 ) 
(m e , m^, m T ), we arrive at the expressions of Mi and \J V The relations 



7e,l = 7 M ,2 ± 90° 



7r,3 



(46) 



V l V V 



must be required in the phase matrices P l v , while Q l are arbitrary. 

The mixing angles 12 , ^23 > ^13 an d the CP phases 5, p, a can be obtained from V 

p}u}qu v p v . 

The numerical results of the parameter space and phenomenological predictions in ansatz (A) 
are shown in Fig. 3. Only 3 free parameters, m 1 , a and j3, are taken as inputs. The experimental 
constraints are the same as those in the general case. The constraint on m l in this ansatz is much 
stronger than that in the general case. One can get m 1 ~ (0.002 — 0.003) eV in Fig. 3. Although 
the number of free parameters has decreased to 3, the numerical results of the mixing angles # 12 , 
#23 and 9 13 still fit the experimental constraints very well. Among them, 6 12 and l3 are still nearly 
evenly distributed in the parameter space, and 9 23 has a very slight preference for being larger 
than 45° . The CP-violating parameters are constrained more stringently. The allowed region of the 
(a,/3) parameters is much smaller: \a\ ~ (45° — 90°) and |/3| ~ (120° — 180°). \J\ can maximally 
reach 0.02 at 9 13 = 8° and 0.03 at 9 13 = 9.6°. The relation p ~ a ± 90° is a good approximation. 
The ratio {fn)ee/{ m )e 1S m ore likely to get a small value than that in the general case. It is only 
allowed in the range (0.002-0.04). Taking (m) e ~ 10" 2 eV, we obtain (m) ee ~ (0.2-4) x 10~ 4 eV. 
This is far beyond the sensitivity of the future experiments. 



5.2 Ansatz (B) 

In ansatz (B), the requirements |CJ = \B[\ and \C V \ = \D V \ are imposed. This ansatz was first 
proposed in Ref. |8j. It is motivated by the mass hierarchy of the charged leptons and the 
experimental fact that the mixing angle 9 23 in the MNSP matrix is about 45°. The relation 
|C,| = \B,\ will lead to \C,\ ~ \mu\i which is compatible with the fact that charged leptons have a 
large mass hierarchy. And the requirement \C U \ = |-D„| can lead to #23 = 45° easily. A detailed 
interpretation for this ansatz can be found therein. Here we reanalyze it by using the latest 
experimental data. 

The solutions for diagonalizing M l and M v in terms of the mass eigenvalues and their ratios 
have been give in Ref. [Sj. We use them for our numerical calculation and show the relevant results 
in Fig. 4. The same inputs and constraints in ansatz (A) are applied to this ansatz. The lightest 
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neutrino mass m 1 is given by m 1 ~ (0.004 — 0.008) eV, bigger than that in ansatz (A). For the 
mixing angles, # 13 > 8.8° and 6 12 < 33.8° hold, and 9 23 is easier to gain a value smaller than 45°. 
As shown in Fig. 4, two thirds of the (^121^23) parameter space is excluded. The constraint on 
the (a, f3) parameter space is still loose and the relation p ~ a ± 90° is also valid. \J\ in this 
ansatz can maximally reach 0.02 at 9 13 = 9.6°, smaller than the maximal value in ansatz (A). 
The prediction for the effective mass of the 0v2/3 decay is totally different from that in ansatz 

(A) . It gives {m)ee/( m )e — 0.1. Since (m) e ~ 0.01 eV also holds in this ansatz, we arrive at 
(m) ee ~ 0.001 eV. We can compare the new results with the old ones presented in Ref. [8]. Since 
the mixing parameters are measured more precisely, most part of the parameter space is excluded. 
Ansatz (B) now is not so favored as before. 

In this section, we have analyzed the UTZT in two ansatze. They have two main different 
features distinguishing themselves from each other. One is the difference of the parameter space 
of the mixing angles. Ansatz (A) is favored in the full (^12^23) parameter space, while ansatz (B) 
is just partly favored. This feature makes ansatz (A), which is a natural assumption of our model 
in section 2, more interesting than ansatz (B). The other feature is the prediction for {m) ee . The 
value of {mn) ee in ansatz (B) is much larger than that in ansatz (A), although both are below the 
sensitivity of the near future experiments. 

6 Conclusion 

The GISM gives vanishing (1,1) and (1,3) submatrices of the 9x9 neutrino mass matrix Ai. This 
is similar to the UTZT which gives vanishing (1,1) and (1,3) elements of the 3x3 mass matrices 
M l v . We have pointed out their similarity and considered their several aspects. The main points 
are listed in the following. 

(1) We have proposed a model based on the discrete Abelian group Z G x Z G to realize both the 
GISM and the UTZT. We reiterate that besides Z§ there are many discrete Abelian groups whose 
direct products can realize both of them. 

(2) We have calculated the UTZT in the general case. Only the normal hierarchy of the neutrino 
masses is allowed by this texture. We obtain the lightest neutrino mass m 1 ~ (0.001 — 0.015) eV. 
The Jarlskog invariant J can maximally reach 0.04 in view of the new experimental results of # 13 . 
The effective mass (m) ee in the 0^2/3 decay can maximally reach 0.001 eV. 

(3) We have compared two ansatze of the UTZT. Ansatz (A) is a natural approximation of 
our model built in section 2, and ansatz (B) is a special case which has been considered in Ref. 
[8]. The mixing angles in ansatz (A) fit the experimental constraints quite well, while in ansatz 

(B) , 8 13 > 8.8° and 8 12 < 33.8° are allowed, and 6*23 < 45° is preferred. Ansatz (B) predicts the 
effective mass (m) ee ~ 0.001 eV in the 0^2/3 decay experiments, while ansatz (A) can only predict 
{m) ee one or two orders of magnitude smaller than that in ansatz (B). 

Finally, we stress that the GISM can avoid the hierarchy problem and is testable in collider 
experiments, and the UTZT agrees very well with current neutrino oscillation data. Both the 
GISM and UTZT can be realized from the same Abelian symmetry due to their similar structures, 
although their uniqueness connot easily be verified in the bottom-up approach of model build- 
ing. Except for the above discussions, there are some other interesting aspects of the GISM and 
UTZT in neutrino phenomenology. One is to discuss possible collider signatures of the TeV-scale 
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right-handed or additional gauge-singlet neutrinos in the GISM, which could be explored by the 
Large Hadron Collider. Another aspect is related to the baryogenesis via leptogenesis, so as to 
account for the cosmological matter-antimatter asymmetry. The two-zero textures of the Yukawa 
coupling matrices and the uncertainty of the scales of M R and may affect how the leptogenesis 
mechanism works in the early Universe. A detailed analysis of these aspects will be done elsewhere. 
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A Simplification of the neutrino mass matrix in the GISM 

A.l General analysis 

The neutrino mass matrix in the GISM is described by a 9 x 9 matrix M given in Eq. (2), where 
M D , Mg, M R and are 3x3 complex submatrices. For physical conditions, one can naturally 
assume that the scale of M D is the electroweak scale and the scale of M s is several orders larger 
than that of M D . To some extent, the scales of M R and are more arbitrary. They can be 
either very high or very small due to different mechanisms. Large mass scales can be regarded 
as the breaking of a certain symmetry at a very high energy scale, similar to the Majorana mass 
matrix of the right-handed neurinos in the type-I seesaw model. And small mass scales may be 
generated from higher dimensional operator after integrating out some unknown heavy fields [25] , 
Small mass scales are also consistent with the 't Hooft's naturalness criterion [26], because the 
conservation of the lepton number is recovered when M R and reduce to zeros. 

Since different scales of M R and may lead to different phenomenological consequences, it 
is necessary to do a general analysis of how M can be simplified in different cases. We denote 

Aft=(M D 0) and M k = (47) 

Obviously, the scale of Mp is several orders smaller than that of M g , and M R yields the masses of 
the right-handed and additional gauge-singlet neutrinos. One can obtain the mass matrix of light 
left-handed neutrinos through a seesaw-like formula 

M v k, -M^M r _1 Md T = -M D (M R - M s M~ l M^)~ l . (48) 

The mass formula in Eq. (48) is the main result in the GISM. It can be further simplified in some 
special cases. However, since Eq. (48) is only valid for M D <C M R — MgM" 1 M§ ', the exception 
should also be considered especially. 

A. 2 Special cases 

For different mass scales of M R and M„, the expression of M v in Eq. (48) can be simplified. For 
the sake of convenience in the following dicussions, we denote M R to be block-diagonalized by a 
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6x6 unitary matrix W as 



W^M^W* = m , (49) 



in which M m and M h are 3x3 matrices standing for the medium and heaviest neutrino masses, 
respectively. Here we consider three typical cases to simplify the mass matrices M u and M R . 
Case (A): s$ M R « M s and < <C M s . Eq. (48) can be simplified to [27] 



M v « MvM^M.M^ME , (50) 



and M R is simplified to 



« M ? Ms . (5!) 

This case has been discussed in Ref. [28J. Since and M R are much smaller than M s , the right- 
handed and additional gauge-singlet neutrinos have nearly degenerate masses and are combined to 
form the pseudo-Dirac particles. Their masses can be not huge and may be testable by the collider. 
For instance, assuming M u ~ 0.1 eV, M ~ 1 keV and M D ~ 10 GeV, we obtain M s ~ 1 TeV. 
Another aspect of this case is the non-unitary effects. Such effects in the mixing matrix are 
approximate to M^M^ 1 . Experimental data show that they are smaller than 0(1) [29]. Due to 
present accuracies for measuring mixing angles, we do not have to consider the non-unitary effects 
in the mixing matrix. We will ignore them in the main body of this paper. 
Case (B): < M s < M R . M R can be simplified to 

M m « M^-MgM^Ms, 

M h « M R . (52) 

This case accommodates a large range of the masses of sterile neutrinos and provides a possibility 
for low scale leptogenesis [3]. One can further discuss the case (Bl): M R < M S M~ X M^ and case 
(B2): M R > M$M~ X M%. In case (Bl), M v can be simplified to Eq. (50); while in case (B2), M v 
can be simplified to [3J 

M v ps -M^M^Ml . (53) 

To derive the tiny left-handed neutrino masses in case (B2), the scale of M R should in general be 
very high, which is similar to the type-I seesaw model. 

Case (C): M R » M s and M M 3> M s . In this case we obtain Eq. (53) and 

M m « max(M R ,M M ), 

M h « min(M R ,M M ). (54) 

The choice of the mass scale of is a little arbitrary except for 3> M s . There is only small 
mixing between the right-handed and additional gauge-singlet neutrinos. 

Since Eqs. (50) and (53) are the typical formulas of left-handed neutrino mass matrices in 
the OISM and type-I seesaw model, respectively, these two models can be regarded as two special 
cases of the GISM to some extent. 
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A. 3 Exception 

Note that Eq. (48) does not hold for M R — M s M~ 1 Mg 1 < M D . This except ion should be considered 
in particular. It can be further divided into two cases: case (D), M R < M D and MgM^M^ 
< M D ; and case (E), M R > M D and MgM^M? > M D , but there is a cancelation that leads to 

Case (D). Since M g is several orders higher than M D , we are led to M R <C M g <C M^. M 
can be simplified by a congruent transformation with a 9 x 9 unitary matrix W. One can write 
out W and W ] MW* as 
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(55) 



respectively. Finally, we obtain M v « M D , which is too heavy to be the left-handed neutrino mass 
matrix. Thus, this case is not interesting. 
Case (E). W and W^MW* are given by 



(\ 






1 





MjM.7 1 



/ o 



yMg'M R " 1 Mp 



R 



1 

M D M R 





(56) 



respectively. One can further derive that the left-handed and additional gauge-singlet neutrinos 
have nearly degenerate masses M v ~ M D M R 1 M g and form the pseudo-Dirac particles. However, 
since M s 3> M D , one has to require that the scale of M R in the GISM be even higher than that in 
the type-I seesaw model, which is unnatural. 
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Figure 1: The comparison of the values between 9 1 and 6 e (left) and that of the values between 9 2 



and 9 T (right). The free parameters 



'3' "m 



are given by Am^ ~ (7.4 
9 12 ~ (33° - 35°) and 13 c 



7.8) x 10~ 5 eV 2 , Am 2 3l 
(8.0° - 9.6°). 



5^, 5 2 , a, (3 and m 1 are used as inputs. The constraints 
(2.4 - 2.7) x 10^ 3 eV 2 , 6 23 ~ (42° - 49°), 



19 





20 




10 


> 




CD 




CO 

1 


5 


O 









CvJ 
CD 



8 8.4 8.8 9.2 9.6 

1 3 [°] 



48 
46 
44 



42 































,' ***** 




. If •■' 










































•- •* . ■ . 







33 



34 

o 12 n 



35 



CO. 



180 



-180 
-180 




90 



180 



-90 



-90 



.... I 



, 7:-.-v.v- 




p[°] 



90 



0.04 



^ 



-0.04 



8 8.4 




(m) e [eV] 



0.02 



Figure 2: The parameter space and phenomenological predictions in the general case. The inputs 
and constraints are the same as in Fig. 1. 
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Figure 3: The parameter space and phenomenological predictions of ansatz (A). Only three free 
parameters a, j3 and m 1 are adjustable. The constraints are the same as in Fig. 1. 
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Figure 4: The parameter space and phenomenological predictions of ansatz (B). The inputs and 
constraints are the same as in Fig. 3. 



22 



